Exact Density of States for finite Gaussian Random Matrix Ensembles via 

Supersymmetry 
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We calculate the exact density of states (DOS) for the three classical and two non-classical Random 
Matrix Ensembles for finite matrix size N using supersymmetric integrals. The 1/N— Expansion 
yields already in lowest order good approximations to the exact result even for small values of N ~ 5. 
We conjecture a connection between the JV— dependence of the oscillating part of the DOS and the 
short-distance behavior of the two-level correlation function. 
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I. INTRODUCTION 

Random Matrix Theory (RMT) has become a well- 
established tool to describe the statistical properties of 
the energy levels of quantum many-particle systems in 
the ergodic regimen. It is based on the assumption that 
(apart from global symmetry conditions as spin rotation- 
or time reversal invariance) the Hamiltonian H of a sys- 
tem with N states is described by a probability distribu- 
tion 



(Hi: 



0. 



(HijH kl ) 



gSikSji- 



RMT derives the statistical properties of the eigen- 
values of H from this distribution of its matrix ele- 
ments. Whereas most of the applications are still cen- 
tered around the three classical Wigner-Dyson (WD) En- 
sembles, the Gaussian unitary, orthogonal and symplectic 
ensembles (GUE, GOE and GSE), there is also interest 
in certain generalizations of the standard cases, leading 
to seven additional symmetry classes, which differ from 
each other by their behavior under a set of discrete space- 
time transformations. Three of them are the chiral ana- 
logues to the WD-Ensembles and four are relevant for 
mesoscopic normal-superconducting hybrid systems!! 

RMT is a soluble description of the statistics of disor- 
dered and chaotic systems in the sense that all n— level 
correlation functions can be computed in principle ex- 
actly. The most important of these correlation functions 
is the two-level correlation function as it contains already 
the main information about the repulsion of neighboring 
levels and discerns the three WD-Ensembles unambigu- 
ously. Several mathematical tools are available to com- 
pute these correlation, functions, e.g. the method of or- 
thogonal polynomialsu, or the sypersymmetric methoda. 
The last method maps the matrix ensemble to a zero- 
dimensional non-linear sigma-model in the limit of large 
matrix size N. Because then the limit N — > oo is implicit 
in all calculations, it is not possible to study directly 
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finite-size effects, i.e. the dependence of the average level 
spacing A(E) at energy E on N . 

Instead one is forced to fix a certain scaling of the 
overall normalization of H with N . In this way the two 
energy scales, the bandwidth 2E max and the mean level 
spacing A(E), are entangled with the matrix size. Either 
the bandwidth is kept at a fixed finite value (macroscopic 
scaling) or the mean level spacing A(E) ~ 1 at a certain 
energy E (microscopic scaling). The first scaling is used 
for the DOS (and yields the semi-circle for N — > oo), the 
second is appropriate to extract the universal level repul- 
sion. The non-classical Ensembles (we will treat class D 
and Class C in this article) show characteristic features of 
the DOS, deviating from the classical ensembles, but 
was calculated only in the microscopic scaling limit. I 
The DOS varies strongly on the scale of the level spacing 
for E ~ 0, and becomes therefore dependent on the posi- 
tion within the spectrum. In this case the usual unfolding 
procedure which normalizes the DOS to a constant over 
many consecutive levels is inapplicable. The question 
arises, whether the new features of the DOS survive the 
macroscopic scaling A(E) ~ 1/N leading to a finite band 
width, which is the physical relevant case for solid state 
applications. It seems therefore advisable to compute the 
DOS exactly for finite N and to study its behavior in the 
limit N — > oo using this formula. The only results for 
finite N were so far obtained via the orthogonal poly- 
nomial method. We will show in the following that the 
supersymmetric method can be used as well while being 
technically simpler and more versatile. 

The outline of the paper is as follows: In section II we 
study the simplest case, the GUE, and explain in detail 
how the exact DOS for finite N can be obtained. We then 
use a saddle-point method to derive a 1/iV-expansion. 
This expansion is not based on a supersymmetric saddle- 
point manifold but to the contrary is effected by exact 
integration over the fermionic variables, which inevitably 
breaks explicit supersymmetry. Taking into account the 
quadratic fluctuations of the remaining bosonic variables 
yields then already excellent approximations to the ex- 
act result for small N and everywhere in the spectrum 
except at the band edge, where the 1/N— expansion di- 
verges. In Section III we give the exact DOS for finite 
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N of the orthogonal and symplectic ensembles as well 
as class C and D. These ensembles are especially simple 
among the non-classical ones because the joint eigenvalue 
distribution can be associated with free fermions on a line 
in an external potentialtl. Nevertheless the calculations 
of the DOS become only simple if the confining part of 
the potential is neglected, corresponding to microscopic 
scaling, thereby loosing all information about the relation 
between A(E) and N which is just given by the confine- 
ment. Therefore we treat Class C/D here together with 
the WD-Ensembles. Regarding the level repulsion they 
belong to the unitary universality class, i.e. without time 
reversal invariance (see below). Our method, of course, 
extends to all other Gaussian Ensembles as well. Sec- 
tion contains the 1/N— expansion around the bosonic 
saddle-points for GOE, GSE and Class C/D. We observe 
a relation between the TV— dependence of the oscillatory 
part of the correction terms and the short distance be- 
havior of the two-level correlation function. Then we 
determine the average level spacing in the band center, 
where the DOS of class C and D deviates strongly from 
the (constant) behavior known from the WD-Ensembles, 
in a way independent from any scaling assumption. Sec- 
tion |v| contains a summary of our results. 

II. THE UNITARY ENSEMBLE 

The ensemble consists of complex hermitian N x N 
Matrices H = H* . We represent the DOS at complex 

I 



(n(z)} 



If we now set z = E + ie with Eel and perform the in-, 
tegration over p and q, we get the exact finite— N results 

e -B 2 

(n(E)) = 2jVjv!V - {H 2 N {E) - H N+1 (E)H N ^(E)) , 

(7) 

where H n (E) denotes the n— th Hermite polynomial. To 
perform an expansion around the limit N — oo, we intro- 

I 



energy z (Imz > 0) as the following expectation value 
(summation convention is understood): 

( n ( z )) = ^ Re J v & ^> H ] <^ Sh ■ (!) 

Here <j> = ipii V^) 7 "; <f> — ( x ifi ) i) ls a vector with N 
bosonic and N fermionic components. We set the vari- 
ance of H to g — 1/2. The action reads 

S H = i4 (Hij - zS^j + trH 2 . (2) 
Integration over H yields an effective action 

S 4 - -iz4<fn + \tvg(4>i4) 2 . (3) 

A Hubbard-Stratonovich transformation to Q— variables 
(the q,p are bosonic and the •& are fermionic), 

«>-(!*)■ (4) 

gives an action bilinear in <p: 

5 Q =i4(Q-z)0 ! +trgQ 2 . (5) 

Now we integrate over the <p— fields and the Grassmann- 
variables fl, i9 exactly to get an action, which depends 
only on the two real variables q and p. The DOS reads 



I 

duce a rescaled energy variable x = E/y/N. With this 
scaling the band width becomes finite and the DOS a con- 
tinuous normalizable distribution in the limit N — * oo. 
After the rescaling (which is done for p and q as well) we 
get for the DOS (after integration over $, d but before 
integration over p and q) 



?e -trgQ 



detg(Q — z) 



N 



ttN 



■ Im 



dpdq 



-qe 



(ip - z) N 
{q-z) N 



N 



2(q-z)(ip-z) 



(6) 



[ [ )) 7T J n/N qe (q~x) N V 2{q-x)(ip-x)J 



which can be written as 

(n(a;)) = / dpdqe- NS(p ' q ' x) f(p, q, x) 



The large prefactor N in the exponent allows for a saddle- 
point approximation of the p, q— integrals in the usual 

(9) 
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way. We get four solutions of the saddle-point equations 

(10) 



q± = x/2±i/2\T2- 
The stability matrix is 
8 2 S 



P± 



-iq± . 



2dq 2 
d 2 S 
2dp 2 
d 2 S 
dpdq 



= l-q\= ±i \J2 - x 2 q± = : Aj 



= 0. 



(11) 
(12) 

(13) 



Figures [j] and § show the location of the saddle-points in 
the complex plane along with their stable directions and 
possible integration paths. The integral (0) converges 




FIG. 1: The saddle-points of q and possible integration paths. 
Solid lines denote the stable directions and dotted lines the 
unstable directions. 
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Re 







contribution surviving in the N — > oo— limit. The action 
S(q-,p+) does not vanish but is purely imaginary 



? (2) 
°SP 

and yields a contribution ~ e lyo sp of modulus 1. This 
contribution disappears in the N = oo— limit, because 
the factor (1 - 2 ( q -x)(i P -x) ) vanis hes for q = q-,p = p+. 
It yields, however, an important oscillatory contribution 
oc 1/N, whereas the (non-oscillatory) corrections to the 
first saddle-point start with order 1/N 2 . The 

final result is 



i(x)) = 



V^X 2 _ (-l^CQS [NSq(x)} 

tt nV2N(2-x 2 ) 



° [ N~ 2 



(15) 



with 

Sq(x) = S 



sp( x ) 



z\Jl - x 2 + 2 arcsin(x/-\/2) (16) 



and 



S' (x) = 2y/2 - x 2 = 27rn (x) 



(17) 



where the zeroth order approximation to the DOS (the 
semi-circle) is written as no(x) = (l/7r)\/2 — x 2 . The 
exact density of states and the 1/N— approximation is 
shown in figure g. From (H) and (0) 

we note that the 




FIG. 2: The saddle-points for p. 

only for Imz > 0. This means, that the integration- 
path of q must pass below the singularity at q = x. 
This forces us to use only the saddle-point q- of the 
g-integration: A path of integration which also crosses 
q+ leads through the region between the two saddle- 
points, where the integrand of @ is large. Therefore 
the correct choice for the path of integration over q is the 
solid line in figure [j]. In contrast it is possible to use an 
integration-path which goes through both saddle-points 
for the p— integral, which is shown in figure [|. 
After the change of variables 



6q 



y/NAZ 



p^p± 



dp 



(14) 



we may perform an expansion in Sq and 5p, which defines 



the 1/N— expansion. The action S 



(i) 

SP 



S(q_,p_) van- 



ishes and yields the well-known semi-circle law as the only 



FIG. 3: Exact DOS of the GUE for N = 7. The dashed line 
is the 1/N— approximation. 

number of maxima of the oscillatory 1/N— term equals N 
and the local width between two maxima scales as the in- 
verse of the DOS in zeroth approximation. The positions 
of the maxima give therefore the locations of eigenvalues 
for a "typical" realization of a random matrix from the 
GUE with one state per maximum on average. Identi- 
fying the distance between to adjacent maxima with the 
average level spacing at x = E/s/N, we get for large N 
the following expression for A(E): 



A(E) 



1 



Nn {E/VN) 



(18) 



which is valid for E not too close to the band edge. The 
l/A^-expansion diverges obviously at the band edge, be- 
cause the saddle-points coalesce and as a consequence the 
stability matrix vanishes. 
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III. 



THE OTHER ENSEMBLES 



Now we introduce the notation 



The elements of the orthogonal, symplectic and Class 
C/D ensembles are conveniently represented as 2Nx2N- 
matrices having a 2 x 2-block structure: 



H = 



A B 
C D 



(19) 



where the N x iV— matrices A, B, C, D fulfill a set of con- 
ditions defining the different ensembles. We set the vari- 
ance of H to g = 1 in the following, which yields a band- 
width of 4 in all cases. 



The orthogonal Ensemble 

The elements of the GOE are real symmetric 2Nx2N- 
matrices: H = H T = H* . To simplify calculations, it is 
convenient to use instead of H the unitary equivalent ma- 
trix M = UHW with U = exp(i(7r/4)<7i)®ljv, where o\ 
is a Pauli matrix acting on the 2 x 2— blocks. If H is real 
symmetric, M fulfills the relation M = a\M T a\. This 
entails the following relation among the block matrices 



M 



A B 

St A T 



(20) 



where A is a hermitian N x N— matrix and B is complex 
symmetric. As in section II, we write for the DOS 

(«(*)> = -L^B e Jv[4>, 0t, M] (x* liXli )e- s , (21) 

with (f> = (xi,X2,ipi,ip2) T , where the £1/2,^1/2 are 
Af— dimensional complex bosonic, resp. fermionic vectors 
and S reads 



S = i0+(M - z)cj)+^trM 2 . 



(22) 



Xij = {x lj ,xl j ,ip lj ,ip 2j ) T , \2j = {x 2jl x* lj ,7p 2j ^ lj ) T 

(23) 

(j is the index of the N— dimensional tensor component) 
together with the transposition 

X*j = (xij,X2j,i>tj,-i>2j) , X% = (^j.^ij.^j. -ipij) ■ 

(24) 

After Integration over M, we get the quartic action in 



S4 = -iXiiXH + ^trg(xiiXH + X2iX2if 
After transformation to the Q-matrix 



(25) 



Q 



/ qi q 2 -&1 -#2\ 

92 qi $2 

i?i $2 ip 

V# 2 #1 ip ) 



(26) 



we obtain 



(n(z)) = ^lmjv[Q} 



qie~ 



detg(Q - z) N 



(27) 



with 



detgQ 

Q 



(ip)' 



■ exp 



.^ 101 - ,262-^62)- |p 



(28) 



and Q = q\- q 2 ql, 9i = 0n?l + #2<?2, 6 2 = 1?i1? 2 - 
Note that the Grassmann-integration contains now quar- 
tic terms in 1?, 1?, which can nevertheless be done exactly. 
The result is (q = qi, r = q 2 q 2 ): 



(n(z)) = -1 Im /f^dr ( 1 - , , + „ „ ^ = ^ J . (29) 

VW/ nN J 7T [(g-z) 2 -rf V [(? - ^) 2 - r] (ip - z) 2 [(q - zf - r] (ip - zf ) K > 

Setting now as above z = E + ie and integrating over p, q and r, we obtain the exact DOS in terms of Hermite- 
polynomials and the error-function: 

, fJP K e~ E ' / H 2N -i(E) 2 -H 2N (E)H 2N _ 2 (E) H 2N -i(E) 2fc! HF* E \\ 

{n{E)) = wnj* { w^)<- liv-ryr (£ m^ {E) - V r 2 erf 71 J J 

(30) 

I 



The symplectic Ensemble with A hermitian and B complex antisymmetric. Pro- 

ceeding as in the case of the GOE (the definition of x, X* 



Here we have 



■ ( ,) at) (31) 



5 



is the same), we find the quartic action 

S 4 = gtrg(xiiX t li + S3X 2l X2^3) 2 , (32) 

where £ 3 = a 3 ®t]y. The corresponding Q— matrix reads 
now 



Q = 

and the DOS reads 
(«(*)) 










02 \ 





q 


-02 




01 


#2 




1P2 


^2 


# x 


ip2 





itrgQ 2 



^Im|p[Q] de f g(Q _ z)jv 



(33) 



(34) 



with 



detgQ" 1 



P 



cxp ( — (ipi6i - ip 2 6 2 



and P = (ipi) 2 - (ip 2 )(ip* 2 ), &i = 0i0i +0 2 2 and 9 2 = 
2 0i. Setting p = Pi and r = P2P 2 , eq,(|3^) reads after 
integration over the Grassmann variables, 



(n(z)) = — — 1m / ar qc 



Uip — z) 2 + r 



N 



(q-z) 



2N 



2N(ip - z) 



(27V + 1)N 



[(ip - z) 2 + r)(q- z) 2 [{ip - z) 2 + r] (q - z) 2 , 



(36) 



and 



HE)) = 



(27V)! 



4 AT 



TV! 



7V-1 



\H 2N (E) 2 - H 2 n+i{E)H 2N - 1 {E)) - —H 2N (E) ^ 



fc=0 



g 2fc (^) 
4 fe fc! 



(37) 



Class D 



with 



The elements of class D are parameterized as 

H = , (38) 

with A hermitian and B complex antisymmetric. The 
quartic action reads (in the notation above) 



1 2 
Si = gtrg ((xiiXij + X2 t x\i)^) 



and the Q-matrix 



Q 



(q x q* 2 -1?! 

92 -<7l ??2 01 

0i J?2 ip 

\02 0i -ip/ 



For (n(z)) we get 



-1 f aie"5 tr gQ 2 
(n(z)) = — lm V[Q] 



ttN 



detg(£ 3 Q - z) 



N 



(39) 



(40) 



(41) 



detg£ 3 Q = 

Q ( 2 ( n n , *n ^ 2 ° 2 ® 2 

w cxp r qm (<zi61 " 9202 + 9292) " qw 



(42) 



and Q = q\ +32^2, ©i = 1M1 - ^2^2, 6 2 = 0102- The 
integrand of contains a pole depending on 9 = gi 
and r — q 2 q 2 . To circumvent it, we have to deform the 
integration path into the complex plane: 



e- i7r ^q - 
e i7r/4 ^ 



ee i7r/4 tanh(g/e) (43) 
-ee- i7r/4 tanh(VF/e) (44) 



for Rez > and e must be small enough: < e < 
(lm z + Re z)/y/2. We find then 
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»(*)> = ^ Im /^dr ^ ^ *J- (l ^-T) . 

AZ)) nN^J n drq ° [( g _ z )2 + r ]W y 1 2[{ q - Z y + r]{ip-zf ] 



and 



-(-l) N e- E ((-1) N (N-1V (Z^(-l) k kl 1 



(46) 



Class C 



We have 



The parameterization reads 

A B 



H 



St -A 1 



(47) 



with A hermitian and B complex symmetric. The quartic 
action is: 



#4 = trg(xi l Xr i S 3 + S 3 X2 l X2 l ) 2 ■ 
with the Q— matrix 

/ q -§ x 2 \ 

—q -0 2 tf x 

i?i i? 2 ipi 1P2 

\i?2 i?i 1P2 -ipi/ 



(48) 



(49) 



-1 /" aie"5 tr eQ 2 
(n(z)) = — Im /2?[Q] iJ 



detg(£ 3 Q 



with 



detg(S 3 Q)- 1 



Pq 2 



(50) 



P f 2 (■ n j.- n • *n a 26262 , ...... 

^■exp ( — (lpiBi +ip 2 B 2 - ip 2 c ) 2) 5I2 - ) ( 51 ) 



andP = (ipi) 2 + (ii5 2 )(i^), 6i = - and 9 2 
??ii?2. With p = pi, r = p 2 p2 the result reads 



W " irN J ir q {q-zf N \ 2[(ip-z) 2 -r]( 



/ \ i ■■' \ 1 i 

(52) 



r] (3 - z) 



{-l) N N\e- E ^ { -ifH 2k {E) H 2N+2 (E) ^ (-l) k H 2k (E) \ 
{n{E))= (2N)l^ [ HME) 2^ fc!4^ + AN ^ ' (53) 

\ k — k — / 



IV. 



SADDLE-POINT APPROXIMATION 



The saddle-point approximation for the four ensembles 
proceeds as in the unitary case in general: Al l four en- 
sembles share the saddle-points q± — x/2 ± iyl — x 2 /4, 
p± = — iq±, tq = 0, (x = E/yN). As in the unitary case, 
the points (q+,p±, r ) have to be avoided by the integra- 
tion path. The fluctuations are quadratic in p and q and 
linear in r around the two remaining SP's (g_,p±,ro). 
Apart from them, we have additional SP's for the GOE 
and the GSE, given in Table |. 

The additional SP can not be reached in the GOE case 
(r SP < 0) but has to be included for the GSE (r SP > 0). 
In class C and D there is an additional saddle-point man- 



ifold at x = 0. However, our results displayed below are 
obtained by confining the analysis to the standard saddle- 
points, which exist for all values of x, thereby assuming 
continuity in the limit x — ► 0. In principle this proce- 
dure could have led to an additional singularity at x = 
similar to the divergence at the band edge. But the com- 
parison with the exact DOS in these cases shows that our 
assumption is indeed correct and the additional saddle- 
points play no role in computing the l/N— expansion. 



We give in the following the results for the l/N- 
expansion of the GOE, GSE and Class C/D. 
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TABLE I: The additional saddle-points for the GOE and the 
GSE 



Ensemble 
GOE 

GSE 


-i- 


q 

X 

2 


r 

4"! 
1-4 


1 



The orthogonal Ensemble 



Lets define the N — oo appro ximation to the DOS 
of the GOE as n (x) = {1/tt)^/1 - x 2 /A. Then the 
1/N— expansion reads to order 1/N 2 



{n (x))=n (x)-±^-^ 



1 



1 



1 



iV 2 1287r 6 n (x) 5 N 2 64TT 6 n {x) 5 
I 



cos [2NSq(x) — arcsin(a;/2)] 



(54) 



with 



S (x) = 2 arcsin(x/2) + xy/l- x 2 /4, (55) 
S' (x) = 27rno(aO, (56) 

analogous to the unitary case. We observe the same fea- 
tures as in the GUE case: The number of maxima of the 
oscillatory part within the band — 2 < x < 2 given by the 
semi-circle is 2N, the total number of levels, so we have 
one level per maximum. The correction diverges close to 
the band edge as expected. Apart from the oscillatory 
contribution there is a non-oscillatory contribution of or- 
der 1/N, However, the information about the level repul- 
sion is encoded in the oscillating part cx cos(2NSq(x)). 
This term is of order 1 /N 2 , which means that the repul- 
sion of levels is weaker in the GOE than in the GUE. 
Formally we can write 

(n(x)) — no(x) + — (non-oscillating term) 

+ (oscillating term) (57) 

with a = 2. For the GUE a = 1. Now the two-level 
correlation function R2 (r) behaves for the Wigner-Dyson 
Ensembles asB: 



i? 2 (r) 



,0 



for r <C 1 , 



(58) 



where r is the distance between two levels on a scale 
corresponding to the average level spacing A(x) ~ 1. 
For the GUE we have = 2 and for the GOE = 1. We 
conjecture therefore 



a 



(59) 



as a relation between the universal parameter 0, which 
characterizes the short distance behavior of the two- 
level correlation function and the exponent a of the fac- 
tor 1/N, which multiplies the oscillating term in the 
1/N— expansion. With this relation we can extract from 
the 1/N— expansion of the one-point function universal 
information about the two-point function, which discerns 
the WD-Ensembles from each other (whereas the DOS is 
the same for all three ensembles in the limit N = 00). 



The symplectic Ensemble 

We have (with the same definition of no (a;) and Sq(x) 
as in the GOE) 



1 



(n(x)) = n (x) + 

1 cos 



1 



N 8ir 2 n {x) 
[NS (x) 4 



sin(x/2)] 



iVs (-l) N V8wino(x)i 



(60) 



Because each eigenvalue appears twice in a 2Nx2N- 
matrix from the symplectic ensemble due to Kramers 
degeneracya, we have only N different eigenvalues and 
therefore N maxima within the band. The oscillatory 
term comes with a 1/V~N prefactor, therefore a = 1/2, 
which gives with (59) the correct universal exponent 
13 = A. 



Class C 



The 1/N— expansion yields 
1 



(n(x)) = n a (x) + 



N8TT 2 n (x) 
1 sin [2NS (x) + arcsin(a;/2)] 



N 



Air 3 xno(x) 



(61) 



Each element of the Class C ensemble has a spectrum 
symmetric with respect to x — 0: To each eigenvalue x 
there is a state with eigenvalue —x. No matrix from Class 
C has an eigenvalue zero. Therefore the DOS vanishes at 
x = 0. The number of maxima in the band equals 2N, 
we have again one (non-degenerate) level per maximum. 
Now the spectral region around x = is different from 
the WD-Ensembles because there is no smooth limit of 
(n(x)} for N — > 00: If one introduces the rescaled variable 
y = 2Nx/tt and considers the microscopic limit x — > 0, 
N ^ co, keeping y finite, the DOS reads 



/ f \\ 1 (■, sm(2iry) \ 



(62) 
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FIG. 4: The GOE and the GSE. N = 10 in both cases. The solid line is the exact DOS, the dashed line is the saddle-point 
approximation. 




FIG. 5: Class C and Class D. N = 5 in both cases. The solid line is the exact DOS, the dashed line is the saddle-point 
approximation. 



which coincides with the result in [||. In terms of y, 
the correction term is independent of N in contrast to 
the GUE, where the oscillatory contribution to the DOS 
vanishes as 1/N for large N. This, however, is not a sig- 
nal of an enhanced level repulsion close to x = in Class 
C. It is merely due to the lack of (approximate) trans- 
lational invariance in the band center, which is caused 
by the mirror-symmetry mentioned above. In the GUE, 
translational invariance is broken only through the term 
proportional to exp(— Y^i=\ ) m the joint probability 
distribution of the eigenvalues A, , which can be neglected 
close to the band center and large N. Therefore the oscil- 
latory structure in (n(x)) is smeared out by the summa- 
tion over all members of the ensemble, whereas in Class C 
the quotient of the probability to find a level at the first 
maximum respective the first minimum of (n) is inde- 
pendent of N. The spectrum is therefore more rigid near 
x = 0. The two-level correlation function #2(2/1,2/1 + y) 
behaves for < y\ , y <C 1 nevertheless as y? and belongs 
therefore to the unitary universality classa. The correc- 
tion term does not modify the macroscopic behavior of 
(n) because the function sin(JVx) / (Nx) tends to zero in 
the Z/2— sense for N — ► 00. That means that the prob- 
ability to find m states in a region of width Ax = 
around zero tends to one for 1 <C to <C N in the limit 



TV — > 00, as in the WD-Ensembles. 

Because the 1/iV-approximation is reliable everywhere 
except at the band edge, we can compute the average 
level spacing exactly in the vicinity of x = without 
recourse to measure x (respective y) in units of the mean 
spacing at a distance of many spacings from zerocl. In 
terms of y (which is exactly related to x and the original 
variable E through N), the average spacing between the 
first and second level > is given by 



(A y (l,2)) = ^(z 4 -z 2 )«l, (63) 



where Zk > denotes the fc— th zero of tan(z) — z. This 
relation is valid for large N but A(j,j + 1) can be com- 
puted via ( |6l| ) for all levels j and j + 1 not too close to 
the band edge and arbitrary N. 
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Class D 

The same mirror-symmetry as in Class C is valid in 
Class D. We have 



(n(x)} = n (x) 



1 



1 



N 8n 2 n (x) 
1 sin[2iV5o(a;) - arcsin(x/2)] 



N 



4:Tr 3 xno(x) 



(64) 



very similar to Class C, but now the DOS is enhanced 
in the band center because there are always two states 
with eigenvalue close to zero. The number of different 
maxima is 2N — 1. Again this feature vanishes in the 
N = oo— limit, because the single additional state at 
x = has measure zero for vanishing average level dis- 
tance A(E) ~ 1/N. As in the Class C ensemble we have 
a = 1 and the universal exponent is (3 — 2. Class D 
belongs therefore with Class C to the unitary universal- 
ity class, which determines the level repulsion even in the 
immediate vicinity of E = e.g. for the second and third 
state away from zero. 

Figures § and || give the exact DOS and the 1/N- 
approximation (dashed lines) for the four Ensembles. 



V. CONCLUSIONS 

We have computed the exact density of states for 
the three Wigner-Dyson Ensembles as well as for Class 
C and Class D by evaluating finite dimensional su- 
persymmetric integrals analytically. In this way the 
N — > oo-limit implicit in most of the previous calcula- 
tions could be avoided. The exact results were then 
employed to test a 1/N— expansion, which differs from 
the usual one because it proceeds only in the bosonic 
sector whereas the fermionic sector is evaluated exactly. 
The 1/N— expansion therefore does not start from a su- 
persymmetric saddle-point or saddle-point manifold. It 
turned out that for all ensembles only a discrete set 



of saddle-points is important, whereas the saddle-point 
manifold at x — appearing in Class C/D is not needed 
for an almost exact computation of the DOS everywhere 
in the spectrum (including x = 0) except the band edge. 
The 1/N— expansion revealed a connection between the 
order a of A -1 multiplying the oscillatory term of the 
DOS and the universal short distance exponent (3 of the 
two-level correlation function, eq.(|59|). In this way we ob- 
tain information about the two-point function from the 
one-point function. We believe that relation (|H^) is al- 
ways fulfilled, if the spectrum is approximately "transla- 
tional invariant", i.e. Rn(xi, X2) depends only on X2 — x\ 
(this is implicit in all calculations using an unfolding pro- 
cedure) . At points in the spectrum where this invariance 
is broken (as at x — for Class C/D) the two-point func- 
tion has to be calculated itself, which can, of course, also 
be done with our method. 

The 1/N— expansion of the DOS yields an unambigu- 
ous determination of the average level spacing A(E) in 
terms of E and N everywhere in the spectrum including 
x = for all ensembles. For Class C/D we conclude that 
the special features of the DOS at the band center vanish 
in the macroscopic N — » 00— limit and are of no relevance 
if the bandwidth is kept finite. 

The other Gaussian ensembles are under current inves- 
tigation. It should also be possible to extend the method 
to non-Gaussian matrix models, see e.g. [Q, which are 
difficult to treat with orthogonal polynomials. Here it 
would be interesting to test the results obtained in the 
large N— limit using singular integral equations by com- 
paring with the exact formulae for finite N and the sys- 
tematic 1/N— expansion. 
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